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We demonstrate how ab initio cluster calculations including the full Coulomb vertex can be done 
in the basis of the localized Wannier orbitals which describe the low-energy density functional (LDA) 
band structure of the infinite crystal, e.g. the transition metal 3d and oxygen 2p orbitals. The spatial 
extend of our 3d Wannier orbitals (orthonormalized Nth order muffin-tin orbitals) is close to that 
found for atomic Hartree-Fock orbitals. We define Ligand orbitals as those linear combinations 
of the O 2p Wannier orbitals which couple to the 3d orbitals for the chosen cluster. The use of 
ligand orbitals allows for a minimal Hilbert space in multiplet ligand-field theory calculations, thus 
reducing the computational costs substantially. The result is a fast and simple ab initio theory, which 
can provide useful information about local properties of correlated insulators. We compare results 
for NiO, MnO and SrTiOs with x-ray absorption, inelastic x-ray scattering, and photoemission 
experiments. The multiplet ligand field theory parameters found by our ab initio method agree 
within ^10% to known experimental values. 
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Many electronic properties of solids can now be de- 
scribed ab initio thanks to the advent of powerful com- 
puters and the development of ingenious methods, such 
as density- functional theory (DFT)^ ^ with local den- 
sity (LDA)^ or generalized gradient (GGA)^ approxima- 
tions, LDA-fHubbard U (LDA+U),^^ quantum chemical 
methods, dynamic mean-field theory,^^ quantum 
Monte-Carlo simulations,^^ and exact diagonalization 
for finite clusters .^^^Ml Nevertheless, for correlated open- 
shell systems with several local orbital and spin degrees of 
freedom, electronic-structure calculations remain a chal- 
lenge. 

Ground-state properties and spectral functions may be 
calculated by exact diagonalization of the many-electron 
Hamiltonian, but this is hindered by the exponential 
growth of the Hilbert-space with the number of corre- 
lated electrons in the system. Exploiting symmetry and 
limiting the number of correlated electronic degrees of 
freedom may enable the treatment of relatively large clus- 
ters, as done in the important case of doped high- Tc 
cuprates, where symmetry in the spin sector allowed Lau 
and coworkers to use clusters with up to 32 Cu02 plack- 
ets, each with a single Cu (1^2 _y2^ and two O p orbitals.!^ 
For local properties, such as excitonic spectra, exact di- 
agonalization for finite clusters becomes much more ap- 
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pealing, as relatively small clusters often suffice. Mag- 
netic anisotropics, ^-tensors, magnetization-dependent 
electron-spin resonance spectra, crystal-field excitations, 
and a manifold of excitonic core-level spectra are usually 
well described using very small clusters. For transition- 
metal and rare-earth compounds, the cluster may often 
be limited to merely a single d- or /-electron cation sur- 
rounded by its nearest neighbor ligands as illustrated in 
Fig.[l] For clusters that small, exact diagonalization is 
equivalent to multiplet ligand- field-theory (MLFT), one 
of the earliest quantum-chemistry methods developed to 
describe the electronic structure of transition-metal and 
rare-earth compounds. MLFT is a highly cost-efficient 
method, able to account for many of the local properties 
and excitonic spectra of correlated materials. 

MLFT calculations traditionally use parameters fit- 
ted to experiments. Despite being a great help for un- 
derstanding and interpreting experimental results, this 
approach is however not completely satisfactory and, 
over the years, numerous theoretical studies have there- 
fore been devoted to obtaining MLFT parameters ab 
initio^^^ Sugano and Shulman^^ calculated the ligand- 
field parameters by constructing single-particle molecular 
orbitals (MOs) as linear combinations of atomic Hartree- 
Fock orbitals and thereby in several cases obtained qual- 
itative agreement with experiments. More often, MO 
theory with a more complete basis is used.^^ After the 
LDA had proven useful not only for s- and p-, but also 
for d- and /- electrons in solids, several authors ob- 
tained MLFT parameters by performing an LDA calcu- 
lation for the cluster and using its Kohn-Sham MOs.'^^^'^ 
Such a calculation breaks the translation invariance of 
the crystal already at the single-particle LDA level, and 
it is necessary to remedy finite-size and surface effects, 
e.g. by embedding the cluster in a set of point-charges 
mimicking the rest of the solid. Such procedures are not 
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FIG. 2: (color online) Left panel: NiO LDA band structure calculated with the large LAPW basis set (thin, black lines) and 
with the small Wannier-orbital basis set consisting of 3 O p (green), 3 Ni d{t2g) (blue), and 2 Ni d{eg) (red) orthonormalized 
A/'MTOs per NiO (thick colored lines). Colors are mixed according to the hybridization between the Bloch sums of the three 
kinds of orbitals. The Fermi level is taken as the zero of energy. Middle panel: Wannier-orbital projected densities of states. 
Right panel: The eight Wannier orbitals. Shown are constant-amplitude contours containing 90% of the orbital charge with 
the color (red or blue) giving the sign. The Ni(0)-centered octahedra have 0(Ni) at their corners. 



well controlled, e.g. depending on the details, the sign of 
crystal- field may change .^^^^ 

Here we use a different route to performing ab ini- 
tio MLFT calculations. Our procedure is similar to 
the method originally devised by Gunnarsson et al*^^^ 
for obtaining the parameters in the Anderson impu- 
rity model and, in the last 15 years, used extensively 
for dynamical-mean-field calculations for realistic solids 
(LDA+DMFT).[SH20l ^^^^^ ^^^^^^ MLFT cal- 

culation by performing a DFT calculation for the proper, 
infinite crystal using a modern DFT code which em- 
ploys an accurate density functional and basis set (e.g. 
LAPWs).^^^ From the (selfconsistent) DFT crystal 
potential we then calculate a set of Wannier func- 
tions suitable as single-particle basis for the MLFT 
calculation-l^^t^ Since the members of such a set are cen- 
tered either on the TM or ligand atoms, we shall call 
them Wannier orbitals. Typically, they are the TM 3d- 
and oxygen 2j9-orbitals which, taken together^ exactly de- 
scribe the DFT 3d- and 2p bands. In general, the set 
should be minimal and span exactly all DFT solutions in 
the energy range relevant for the property to be calcu- 
lated. It is important that this set contains sufficiently 
many ligand orbitals to make the correlated TM orbitals 
well localized, i.e. the TM d-orbitals should not have 
tails on any other atom. This localization allows one to 
restrict the many-electron calculations of local proper- 
ties to a single TM site plus its ligand neighbors. Hence, 
in the current method, there are no embedding errors, 
except those arising from truncating the single-particle 
basis to include only the Wannier orbitals on the cluster. 

In the following we introduce the method by the ex- 
ample of the late transition-metal oxide NiO with config- 



uration 3d^. In Sect. II, we show that similar results can 
be obtained for middle and early transition- metal oxides, 
specifically 3d^ MnO and 3d^ SrTiOs- 



In Sect. HI, we compare with results obtained by 
several different experimental techniques: (A) 2p x- 
ray absorption (XAS), a charge neutral excitation of a 
transition-metal 2p core electron into the 3d shell. (B) 
2p core level x-ray photoemission (XPS) from Ni impuri- 
ties in MgO. (C) Inelastic x-ray scattering (IXS) of core 
to valence excitations, a technique similar to XAS. We 
specifically show 3p core-electron excitations into the 3d 
shell. (D) Inelastic x-ray scattering of d-d excitations. 
The experiments presented for these materials are rela- 
tively well understood, so that the comparison with our 
new ab initio results constitute a critical test of the theo- 
retical method. At the end of the paper, we conclude. In 
appendix [A] we provide information on the different ba- 
sis sets or Wannier orbitals used, as well as the meaning 
of the different occupation numbers and the concept of 
formal valence. Details of the calculations, including nu- 
merical values of several MLFT parameters obtained ab 
initio, may be found in Appendix [B] A discussion of the 
double counting of interactions in the LDA and MLFT 
calculations may be found in Appendix [C] In appendix 
|D]we show how Ligand orbitals can be obtained in gen- 
eral symmetry from the O 2p orbitals, with the use of 
blocktridiagonalization of the orbital basis set. This is 
an essential ingredient which makes these calculation nu- 
merically efficient. Appendix [E] contains a short note on 
the exact diagonalization routines. 
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I. OBTAINING THE MLFT PARAMETERS 
FROM THE LDA BY THE EXAMPLE OF NiO 

In this section we introduce the ab initio MLFT 
method by the example of NiO. We will discuss the differ- 
ent steps taken in order to obtain the MLFT parameters. 
First we discuss the LDA procedure used to obtain the 
potential, Wannier functions and tight binding parame- 
ters. Next we discuss the meaning of the different one 
electron parameters. In the last part of this section we 
discuss many body parameters, i.e. the Slater integrals. 

We start our ab initio calculations with a conventional 
charge-selfconsistent LDA calculation for the experimen- 
tal crystal structure. NiO has the rock-salt structure 
in which each Ni atom is surrounded by six O atoms 
in cubic symmetry, and vice versa. We used the linear 
augmented plane wave (LAPW) method^^ as imple- 
mented in WlEN2K.ESlThe resulting LDA band structure 
is shown along the symmetry lines of the face centered 
cubic (fee) Brillouin zone in Fig. [5] It is not very differ- 
ent from the band structure obtained and discussed forty 
years ago by Mattheiss^^ who used Slater exchange and 
a non-selfconsistent potential construction. The three O 
2p bands extend over 5 eV, from -8.2 to -3.2 eV below 
the Fermi level. The five Ni 3d bands consist of three t2g 
bands extending from -3.0 to -0.9 eV and two eg bands 
extending from -1.4 to -1-1.3 eV. The bottom of the Ni 4s 
band is 1.5 eV above the Fermi level and at the F point. 
As pointed out by Mattheiss, the reason why the 4s band 
is above the 3d bands and thus empty, while it is half- 
full in elemental Ni, is that strong hybridization with the 
O 2p band pushes it up (and the 2p band down) in the 
oxide. 

Within the LDA, NiO is a metal, in strong contrast to 
experiments where NiO is found to be a good insulator 
with a room-temperature resistance of ~ 10^ Qcm and an 
optical band-gap of about 3.0-3.5 eV.I^^H^This is one of 
the most noticeable failures of the LDA. However, for the 
current paper, this is not a problem. Although the LDA 
cannot reproduce the correct electronic structure near 
the nickel atom, the minimal set of localized Ni d and O 
p orbitals which together span the low-energy solutions 
of Schrodinger's equation for the LDA crystal potential 
exactly, i.e. the 5 + 3 = 8 bands in Fig.[2j is expected to 
constitute a good single-particle basis set for many-body 
calculations. 

In order to prevent double counting of the multipole 
part of the Coulomb interaction, we constrain the self- 
consistent LDA potential to be spherically symmetric in- 
side non-overlapping muffin-tin (MT) spheres (see Ap- 
pendices |B] and [C|, but allow it to be general in the MT- 
interstitial; it is a so-called warped MT potential. For 
this potential we generate a basis set of 8 localized TM 
d and O p orbitals per cell which span the 8 bands ex- 
actly. Since these bands do not cross any other bands in 
NiO, this can be done by projection of the LDA LAPW 
Bloch states onto Wannier functions choosing band- and 
k-dependent phases which make the Wannier functions 



atom-centered and localized. For an oxide like SrTiOs, 
the TM d and O p bands do cross and hybridize with 
other bands far away from the Fermi level as can be seen 
in Fig. [5] near avoided crossings it is therefore not clear 
which of the bands to project on to. Moreover, one might 
want to go beyond perfect crystals. Rather then using 
projection, we generate the minimal basis set of localized 
orbitals directly by using the A/'th-order muffin-tin orbital 
(7VMT0) method.l^St^ll This method solves the problem 
exactly by multiple scattering theory at + 1 chosen en- 
ergies, followed by A/'th-order polynomial interpolation in 
the Hilbert space, but only for a superposition of spheri- 
cally symmetric short-ranged potentials (to leading order 
in the potential overlap) . We must therefore first perform 
the overlapping muffin-tin approximation (OMTA)^^ to 
the warped MT potential, i.e. by least-squares minimiza- 
tion determine the radial shapes of the overlapping po- 
tential wells and the common potential zero. 

The resulting basis set of five Ni d plus three O p 
ATMTOPI^ per ceh with the A/" + 1 = 2 energies, 
eo = —5.2 and ei = —1.2 eV, produces the eight col- 
ored, thick bands in Fig.|2j Within the width of the 
line they are indistinguishable from the LAPW bands. 
Hence, the A/'MTO minimal basis set for the OMTA 
to the warped potential is a highly accurate represen- 
tation of the large LAPW basis set for the low-energy 
states, but many times more efficient. Our Wannier or- 
bitals are symmetrically orthonormalized A/'MTOJ^^lEIl 
and the colors indicate the relative O p, Ni d{t2g) , and 
Ni d{eg) Wannier-orbital characters. The middle panel 
of Fig.|2] shows the partial density of states projected 
onto these Wannier orbitals. Compared with the com- 
monly used projection onto truncated partial- waves in- 
side a MT sphere, our projection has the advantage of 
leading to well-defined occupation numbers because it 
is onto a complete, orthonormal basis set of localized, 
smooth orbitals. Our projection also takes care of the 
O p(tiu) character which flows into the neighboring Ni 
MT sphere due to the large size of the Wannier O p or- 
bital. In this regard, it should be remembered that a 
MT sphere is not chosen to give a good representation 
of the charge density, and hence of the occupied Wan- 
nier orbitals, but of the potential. Since the latter has 
an envelope function which for rocksalt-structured NiO 
is essentially the Coulomb potential from equal, but op- 
posite point charges on identical cubic sublattices, Ni 
and O have similar sized MT spheres. This makes it 
necessary for the wave-functions resulting from a MT- 
based method for solving Schrodinger's equation such 
as LAPW, to carry the partial-wave expansions much 
further than to p or typically to / ~ 8, because the 
outer part of the O p orbitals, for instance, are being 
expanded around the Ni sites. Nevertheless, with appro- 
priately normalized partial waves, projection of the den- 
sity of states does give similar results as with Wannier 
pd-orbitals. 

The eight Wannier orbitals, Wi (r), are shown on 
the right-hand side of Fig. [2] as those surfaces where 
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\wi (r)| = const and which incorporate 90% of the charge, 
\wi (r)l^ (fir = 0.9. The red/blue color of a lobe gives 
its sign. As one can see, the Ni d orbitals are extremely 
well localized. This is a necessary condition for several 
many-body models which implicitly assume such an or- 
bital basis set, for example the Hubbard model which 
neglects all off-site Coulomb correlations. In order to 
visualize the localization of the Ni 3d Wannier orbitals 
at a more quantitative level, we computed the effective 
radial wave-functions for the t2g and Cg orbitals by mul- 
tiplying with the corresponding spherical harmonics and 
averaging over all solid angles. These radial functions are 
compared in Fig. [s] with that of a Ni atom in the d^ con- 
figuration as calculated with the Hartree Fock method.!^ 
Although there are slight differences, the agreement is 
astonishing. The local Ni d Wannier orbitals in NiO are 
rather similar to atomic Ni wave- functions. Note that 
the atomic Ni d radial function depends on the filling of 
the (i-shell, but is rather insensitive to the filling of the 
As shell. The atomic radial function shown in Fig. [3] is 
calculated for a Ni^+ [3d^4s^) configuration, but would 
be practically the same for a Ni+ (3d^4s^) or neutral Ni 
(3(i^45^) configuration. 
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FIG. 3: (color online) Angular averaged radial wave functions, 
R{r), for the Ni Cg and Ni t2g Wannier orbitals compared 
with the Hartree-Fock radial wave function for a Ni^^ ion in 
a 3d^ configuration. For Ni+ 3d^s^ and Ni 3d^s^, the radial 
functions are similar. The distance to the nearest oxygen is 
2.09 A, which is consistent with the sum of the ionic radii of 
0.72 A for Ni^+ and 1.40 A for 0^~. The inset shows the Slater 
integrals Eq. ([T]) for the multipole Coulomb interactions. 

Since we have chosen not to include Ni s orbitals in the 
minimal basis, it does not describe the high- lying, empty 
Ni 4s band which has anti-bonding O 2p character. The 
corresponding Ni 4^ bonding character of the O 2p-like 
band is however completely taken care of by including 
(downfolding) the Ni 45 character into the tails O 2p 
Wannier orbitals, as is seen in Fig. [2] In the bottom right- 
hand panel one can see how Ni 4s character is added at 
the tip of each lobe of the O 2p orbital, such that the 
outermost 4s radial lobe expands the tip of the 2p lobe, 
while the remaining inner radial 4s lobes of alternating 
sign cause the 2p lobe to tail-off in an oscillating manner. 
The shape of the O 2p Wannier orbital is of course also 



influenced by the requirement that it be orthogonal to 
the Ni 3d Wannier orbitals. 

The A/'MTO method is particularly useful when a 
real-space tight-binding representation of the Hamilto- 
nian is needed.^'' Both the orthogonal Wan- 

nier functions as well as the corresponding tight-binding 
representation of the Hamiltonian in this basis set 
are directly available in the A^MTO formalism. De- 
tails on the A^MTO method can be found in previous 
publications^^ and Appendix |b] Here we would like to 
stress that the Wannier orbitals used within this paper 
are not constructed by maximally localizing the Wan- 
nier functions,'^ but their extend is a result of symmetry 
constraints. This leads to orbitals that are not always 
maximally localized, especially in the details of the tails 
of these orbitals. The Ni d Wannier orbitals obtained 
by 7VMT0 do show a very large overlap with atomic or- 
bitals, which might well be larger than the overlap one 
might find between atomic orbitals and maximally lo- 
calized Wannier orbitals. It is the agreement between 
our Wannier orbitals and atomic orbitals which makes 
the method successful. An alternative method to obtain 
good Wannier orbitals for correlated model calculations 
could be to maximize the overlap of the Wannier orbital 
with atomic orbitals. 

Although only the Ni d bands in Fig. [2] are partly oc- 
cupied, inclusion of O p orbitals in the basis is important 
for describing spectroscopy. In photoemission, for exam- 
ple, the removal of a TM d electron can lead to a transfer 
of charge from the O to the TM atom. This dynamical 
screening would not be captured on a basis of only TM 
d orbitals. Multiplet Crystal Field Theory (MCFT), i.e. 
local calculations using a basis of only TM d orbitals, 
which have antibonding O p tails fixed to them, can be 
useful in many other cases, for example for calculating 
magnetic anisotropics. In this paper, however, we con- 
centrate on Multiplet Ligand Field Theory (MLFT) and 
explicitly include the O p orbitals in the basis set. 

Until now, we have considered the infinite crystal and 
have calculated the one-electron potential in the local- 
density and warped muffin-tin approximations. For that 
potential we have derived a set of localized O p and Ni 
d Wannier orbitals which exactly describe the O p and 
Ni d bands, as well as the corresponding tight-binding 
Hamiltonian. We now use these orbitals and this Hamil- 
tonian for the NiOe cluster (Fig.jl]). The band structure 
thereby reduces to the O p like tt- and <j-levels and the Ni 
d like TT*- and a* -levels shown in the central part of the 
left-hand panel of Fig.|4j plus some O p levels which do 
not couple to the Ni d levels and are therefore not shown 
in the figure. For comparison, we repeat from Fig. [2] the 
crystalline density of states projected onto the O p, Ni 
d (t2g) and d (cg) orbitals. In the following we discuss the 
formation of these simple cluster levels before we consider 
calculating many-electron multiplets. 

The NMTO Ni d Wannier orbitals have the on-site en- 
ergies Ct^g = —1.55 and e^^ = — 1.05eV with respect to 
the Fermi level. Had the potential been spherically sym- 
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FIG. 4: (color online) Left panel: Orbital energy level diagram for the NiOe cluster on the same energy scale as the LDA 
band-structure for the solid (Fig. [2]) shown in the middle panel. The Fermi level is the zero of energy. Right panel: Constant- 
amplitude contours of the Ni d Wannier orbit als and of the Ni-centered ligand orbit als. The latter are symmetrized linear 
combinations of the O p Wannier orbit als. 
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metric within the range of the Ni d orbit als, the Cg and t2g 
radial functions in Fig. [3] would have been identical, and 
their levels degenerate with energy = ^^t2g + f^e^ = 
— 1.35 eV. The crystal- field splitting, lODq = ~^t2g — 
0.5 eV, is basically due to the fact that Cg and t2g or- 
bitals point respectively towards and between the nearest 
oxygen neighbors, which are negatively charged. (The 
notation l{)Dq for - e^^^ is standard in MLFT.I^ 
In the conventional ionic picture, 2 electrons are trans- 
ferred from each neutral Ni 3(i^4s^ atom to each O atom, 
where they complete the 2p shell. Hence, the crystal-field 
splitting would be the radial matrix element of the non- 
spherical part of the electrostatic Madelung potential 

^ [I40 (f) + 7(5/14) [1^44 (f) + 1^4-4 (f 
±2 charges. However 

shown in Fig. 3 which yields: (r^) ~ ^0.71 , the 

splitting is merely ~ 0.3 eV, i.e. ~ 0.2 eV too small, and 
this is even an overestimate because the charge trans- 
fer from the 4s to the 2p band is reduced by covalency. 
Note in passing, that had we taken the anisotropy of 
the electrostatic potential to be the one caused by the 
LDA charge density and the protons outside the Ni MT 
sphere, we would have gotten merely 0.01 eV. This is so 
because the Ni MT radius of 1.10 A is larger than the 
Ni^+ ionic radius of 0.72 A and thus cuts off part of the 
oxygen charge density (remember: MT spheres are de- 
signed to describe the potential and not the charge den- 
sity). Hence, the anisotropy felt by the different angular 
behaviors of the Cg and t2g orbitals can at most only ac- 
count for half the calculated crystal- field splitting. Next, 
we now turn to the different radial behaviors (Fig.[3|. 
These are mostly due to the requirement that the Ni eg 
and t2g Wannier orbitals be orthogonal to the nearest O 



2p orbitals. The eg radial function changes sign for in- 
creasing r, while the t2g radial function merely decays. 
At short distances (r < 0.7 A), the normalized eg radial 
function is therefore larger than the normalized t2g func- 
tion. Since the maxima of the two radial functions occur 
where the radial potential, '^Ni (^) + is huge and 

negative, the higher eg maximum causes a lower poten- 
tial energy, opposite to what is needed to explain the size 
of the crystal-field splitting. In the end, it turns out that 
the ~0.5 eV crystal-field splitting is a result of not only 
the potential acting on the different angular and radial 
wave-functions, but also due to the kinetic energy. The 
eg orbitals overlap more with the O p orbitals than the 
t2g orbitals, whereby orthogonalization increases the ki- 
netic energy more for the former than for the latter. For 
the calculation of the MLFT parameters it turns out to 
be important to treat all these interactions on an equal 
footing. 

We now continue our explanation of the orbital en- 
ergy level diagram in Fig. [4j this time starting from on- 
site energy = —4.74 eV of the O 2p Wannier orbitals 
marked on the right-hand side. The ionic energy is thus 
^d — ^ 3.4 eV. Since we have chosen not to include Ni 
45 orbitals in the basis set, the bonding Ni s character 
has been downfolded into the O p orbitals so that the 
—4.75 eV includes a downwards shift of about 1 eV from 
Ni s covalency. The NiOe cluster contains 6 x 3 = 18 O 
p orbitals, but not all linear combinations can interact 
with d orbitals on the central Ni site. Hence, the basis 
set for the MLFT calculations can be greatly reduced by 
including only those linear combinations which do cou- 
ple, the so-called ligand (L) orbitals.!^ The reduction of 
the Hilbert space by use of Ligand orbitals is crucial for 
efficient MLFT calculations which is explained in more 
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detail in Appendix |DJ For each TM d orbital there is ex- 
actly one such linear combination. The right-hand panel 
of Fig. [4] shows the 5 Ni d orbitals together with the 
5 Ni-centered L orbitals of the same symmetry. There 
is an important difference between the L t2g and L eg 
orbitals: Whereas the L t2g orbitals are bonding (same 
color) between nearest O p Wannier orbitals, and thus 
give rise to a substantial 0-0 <j-like bond charge, the L 
eg orbitals are anti-bonding (different color). As a result, 
the energies of the L t2g and eg orbi tals a re respectively 
Tpp = ppcr — PP7T below and above €p^^^ 

We finally complete the level diagram by including the 
covalent hopping integrals Vt^g = pdirxA/ \f\ = 2pd7T and 

= -^pda X 4/V4 = -V3pda - Spdii^^ between 
the L p and TM d orbitals of respectively t2g and eg 
symmetry. The t2g hopping gives rise to an L p-\ike tt 
and a TM d-\ike tt* level, and the eg hopping gives rise 
to an L p-like a and a TM d-\ike a* level. It is these 
TT* and cr* levels which in the solid broaden into Ni t2g 
and e^ bands. The a* level is close to the Fermi level 
in the LDA and this indicates that the cr* orbital is half 
full. The TT, cr, and tt* orbitals have considerably lower 
energies and are fully occupied. 

Our MLFT calculations include Coulomb correlations 
beyond the one-electron mean-field potentials discussed 
so-far, but only among the TM d orbitals. Arguments 
for treating the L p orbitals as well as their Coulomb 
repulsion with the TM d orbitals at the LDA level, are 
that the L p orbitals are fairly delocalized and that they 
are almost fully occupied. As an example, we can safely 
neglect correlation in an event where two holes meet on 
a single oxygen atom and scatter. The Coulomb correla- 
tions are responsible for the multiplet structure, and we 
keep them among the Ni d orbitals, but make a distinc- 
tion between the spherical ([/, A) and the non-spherical 
repulsions. 

The spherical part of the Coulomb repulsion, often 
parametrized by [/, is strongly screened in a solid. If a 
Ni d electron is removed, there will be a charge-flow into 
the Ni 45 orbital, for example, which reduces the energy 
cost of such an excitation. Although several calculations 
of the screened U have been presented in the past,^^ 67 77 
we fit U such as to obtain the best agreement between 
our M LFT c alculation and the experimental multiplet 
spectra.l^^^ The parameter A is the orbitally averaged 
(spherical) part of the difference between the on-site en- 
ergies of the Ni d orbitals and the L p orbitals at a filling 
of 8 electrons in the Ni d shell. In the LDA, as well as 
for the ground-state found in our MLFT calculations, the 
Ni d occupation exceeds 8 due to the covalency with the 
oxygens. The relation between A, as defined in MLFT 
calculations, and e^^ — e^, as obtained from the LDA, is 
rather non-trivial and we shall therefore treat, not only 
[/, but also A as an adjustable parameter. In the fore- 
seeable future, it should be possible to calculate U and 
A from first principles. 

The non-spherical parts of the Coulomb interactions 
we can easily calculate because the multipole interactions 



between two d electrons are hardly screened. For exam- 
ple, the Coulomb repulsion between two dx2_y2 electrons 
is obviously larger than that between a dx2_y2 electron 
and a (is^s-i electron, but to screen this difference re- 
quires electrons with high angular momentum around 
the Ni site; a Ni 4^ electron, for instance, could not do 
it. Also electrons on neighboring sites are inefficient in 
screening the multipole because it decays fast (oc r~^~^). 
It has been shown that neglecting any screening of the 
multipole part of the Coulomb interaction gives reason- 
able agreement between theory and experiment.'^ Also 
in the present work, we shall neglect any screening of the 
multipole part of the Coulomb interaction and shall find 
reasonable agreement with experiments. 

Multipole interactions are the cause of the Hunds-rule 
energy. For example, two dx'2-y2 electrons must have dif- 
ferent spins, whereas two electrons in different d orbitals, 
and hence less repulsive, may be in a spin-triplet state, 
as well as in the spin-singlet state. Experimentally it 
has been shown^^ that the multipole interactions of the 
Coulomb interaction, are reasonably well approximated 
by assuming that the d orbitals have the pure-angular- 
momentum form: R{r)Y2ra (r)- The inset in Fig. [s] is a 
table of the values of the Slater integrals obtained using 
the Ni^+ ionic radial function, R (r), as well as the radial 
functions obtained by averaging the Ni t2g and eg Wan- 
nier orbitals over solid angles. The Slater integrals for d 
orbitals are: 

F(''^ = j j ^R\r^)R\r^)rldr,rldr2. (1) 

where r< = min(ri,r2), r> = max(ri,r2), and k = 
0, 2, or 4. The definitions of U and the Hund's rule ex- 
change, Jh-, vary: The average repulsion between two 
d orbitals is: Uav = F^^^ - ^(F^^^ + F^"^^). How- 
ever, in order to discuss the Mott gap, one uses the 
energy difference between the lowest multiplets of dif- 
ferent configurations and that has lead to the definition: 
U = F^^^ + + The Hund's-rule exchange 

can either be defined as: Jh = ^(i^^^^ + F^^"^)^ or as: 
= H^^^^ + ^F^^\ The bare F^^) as calculated 
from the Wannier orbitals is of the order of ~25 eV. This 
is clearly much too large because the monopole part of 
the Coulomb repulsion is strongly screened. The values of 
F^^^ and F^^^ are respectively ~ 11 and ~ 7 eV, in good 
agreement with experimental values, as we shall see. The 
multiplet interactions are quite large and lead to a multi- 
plet splitting of the Ni-d^ configuration of about 7.5 eV, 
which is the energy difference between the ground- 
state configuration and the highest excited singlet of 
character. This is larger than the Ni-d bandwidth and 
therefore not a small energy. 

We will compare our results to several experiments, in- 
cluding core level spectroscopy. Once a core hole is made, 
the interaction between the core and valence electrons 
becomes important. Here again we will make a distinc- 
tion between the multi- and monopole interactions. The 
monopole interactions U2p^M and Usp^sd will, like for the 
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valence states, be taken from fits to experiment l ^^ * ^^ * For 
tiie multipole interactions, i.e. tiie Slater integrals F"^^^^ 
Gp ^ and we again assume that screening can be 

neglected, which allows us to directly calculate these in- 
tegrals from the core and valence Wannier orbitals. The 
core Wannier orbitals are equivalent to atomic wave func- 
tions since they have no inter-site overlap. It is important 
to use a scalar-relativistic method for the calculation, as 
well as to calculate the core wave- functions for the final 
state occupations, which requires an open shell calcula- 
tion. We used an atomic Hartree-Fock code to obtain 
these core wave- functions,'^ but any open-shell, scalar 
relativistic method should give similar results. Specific 
values of the Slater integrals can be found in Table [l] in 
Appendix [C] 

We now have all ingredients needed to perform MLFT 
calculations of experimentally observable quantities. But 
before we do this, we will introduce similar ligand-field 
calculations for MnO and SrTiOs. This allows us to com- 
pare oxides of early, intermediate and late transition met- 
als and show that the method is likely to apply to a range 
of correlated transit ion- metal compounds. 



II. NiO, MnO AND SrTiOs 

In Fig.|5]we shown from top to bottom the LDA band- 
structures of NiO, MnO, and SrTiOs calculated in the 
same way as in Fig. |2] and explained in the previous sec- 
tion, with details given in Appendix [C| Whereas NiO 
and MnO have the fee rocksalt structure, SrTiOs has 
the simple cubic (sc) perovskite structure in which the 
Sr cube has Ti at its body center and O at its face cen- 
ters; in MLFT we treat the TiOe cluster. Going from 
NiO to MnO, the TM-electron and -proton counts are 
both reduced by 3, whereby the d-band filling gets re- 
duced from d^ to d^. Concomitantly, we see that the 3d 
bands move up in energy relatively to the 4s and O 2p 
bands, by about 1.5 eV. The p and e^ bandwidths as 
well as the Qg-t2g splitting decrease, presumably due to 
the increased ionicity, Cd — ep. Going finally to SrTiOs, 
the TM-electron and -proton counts are further reduced 
by 3, but due to the change of stoichiometry, the nomi- 
nal d-band filling is now reduced to d^ rather than to d^. 
SrTiOs is a band insulator and the LDA bandstructure 
shown in the bottom panel agrees with the ionic config- 
uration Sr^+Ti^+(0^~)s : We see nine fuh O 2p bands 
separated by a 2 eV gap from the three empty Ti 3d {t2g) 
bands. The latter are separated by a small gap from the 
two Ti 3d(e^) bands which overlap the two Sr 4d(e^) 
bands and the bottom of the Ti 4s band. The three Sr 
4d {t2g) bands are pushed up in energy by covalent pda 
interaction with the 12 nearest oxygen neighbors^^ and 
thus lie more than 8 eV above the Fermi level. Due to 
the different structure and stoichiometry of SrTiOs, its 
bands are quite different from those of NiO and MnO. 

The agreement between the O p-like and TM d-\ike 
bands obtained with the LAPW method and those ob- 
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FIG. 5: (color online) Same as Fig.[2] but for fee NiO (d^), 
fee MnO (d^), and se SrTiOs (d°). On the right-hand side, 
the oxygen orbitals are not shown. 

tained with the minimal basis set of A^MTOs is almost 
perfect for NiO and MnO. The agreement is also satisfac- 
tory for SrTiOs, although near the bottom of the Ti t2g 
band, and at various places in the O 2p band, small dis- 
crepancies may be detected. These are most likely due 
to the OMTA causing a slightly inaccurate description 
of the hybridization with the high-lying Sr Ad and Ti 45 
bands. 

From the right-hand side of the figure, we see that for 
all three materials the TM 3d Wannier orbitals are very 
well localized. This is a necessary condition for using 
them in MLFT. We do not show the O 2p orbitals as in 
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FIG. 6: (color online) Angular averaged TM t2g Wannier- 
orbitals in NiO, MnO, and SrTiOs. The distance to oxygen is 
2.09 in NiO, 2.21 A in MnO, but only 1.95 A in SrTiOs. The 
ionic radii are 0.72 for Ni^+, 0.80 for Mn^+, 0.90 for Ti^+, but 
only 0.68 A for Ti^+. 



Fig.|2l but had we done so for SrTiQq, we would have 
seen not only bonding Ti 4s character of the p orbital 
pointing towards Ti, but also bits of bonding Sr M and 
55 characters of the two other p orbitals, which point to- 
wards Sr.^^ The good localization of the Ti Cg orbitals is 
related to the feature seen in the left-hand panel around 
8 eV, that the A^MTO Ti eg band interpolates smoothly 
across the avoided crossing of the LAPW Ti eg and Sr 
e^ bands. Had this not been the case, the Ti Cg Wan- 
nier orbitals would have been long-ranged. Hence, we 
can construct a minimal set of localized A^MTOs, even 
when bands described by the set are crossed by and hy- 
bridizes with other bands, provided that we can choose 
the + 1 expansion energies outside the range of those 
other bands. For SrTiOs, we used eg = —2.6 and ei = 1.5 
eV. 

The 3d Wannier orbitals for the three oxides are very 
similar; they merely contract along the 3d row of the 
periodic table. This is seen when going from the bot- 
tom to the top in the right-hand panel of Fig.[5| and 
even more clearly in Fig. |6] where we show the angular- 
averaged t2g Wannier orbitals. The well-known reason 
for this orbital contraction is as follows: Upon proceed- 
ing one step along the TM row, a proton and an elec- 
tron are added. The electron incompletely screens the 
attractive potential from the proton seen by another va- 
lence electron, and as a result, the one-electron potential 
becomes deeper and deeper upon proceeding along the 
series, until the d shell is full and the screening is com- 
plete. Since this mechanism is atomic, the shapes of our 
3d Wannier functions are fairly robust and the chemistry 
merely changes tails -and thereby normalizations- a bit. 
This is what we saw in Fig.|3] For that reason, the con- 
traction seen in Fig.|6]of the t2g radial functions -which 
are less influenced by O than the Cg functions- closely 
follows that of the 2+ ionic radii, which are 0.72 (Ni^+), 
0.80 (Mn2+), and 0.90 A (Ti^+j, in the sense that at the 
respective radius, all three radial functions have about 
the same amplitude. This happens although the Wan- 
nier orbitals are calculated for the real solids, which in 



the case of SrTiOs have a Ti-0 distance far smaller than 
expected from the size of the Ti t2g function in Fig. [6] 
That the chemical binding of SrTiOs is different than 
those of NiO and MnO also becomes clear by adding the 
accepted ionic radius of O^" (1.40 A) to the TM^^ radii 
given above, obtaining: 2.12, 2.20, and 2.30 A for the 
TM-0 distance in respectively NiO, MnO, and SrTiOs- 
The real distances in NiO and MnO are nearly the same, 
but much smaller (1.95 A) in SrTiOs- 

The short Ti-0 distance is of course reproduced by 
using the Ti^+ radius of 0.68 A which corresponds to the 
band-structure configuration Sr^+ Ti^+ (0^~)3 = Sr4d° 
Ti3d^ (0 2p^)s- This ionic picture of the binding seems 
to neglect the Ti-0 and Sr-0 covalencies predicted by the 
LDA, i.e. the fact that there is a considerable amount of 
Ti-3d and St-M partial- wave character in the O 2p bands. 
But this is only apparently so: The Ti 3d and Sr Ad radial 
functions are essentially exponentially increasing because 
they solve the respective radial Schrodinger equation for 
O 2p-band energies, which are way below those of the Ti 
3d and Sr 4d bands. Hence, these partial waves simply 
complete the shape of the O 2p Wannier orbitals inside 
the Ti and Sr MT spheres. 



III. COMPARISON WITH EXPERIMENTS 

In order to test the quality of MLFT calculations using 
the LDA Wannier orbitals, we now present a comparison 
between theory and experiment for several established 
spectroscopies, which show TM excitons. Such locally 
bound states are represented well within the small clus- 
ter used in the MLFT. The materials considered, namely, 
NiO, MnO and SrTiOs, are insulators, thus justifying the 
theoretical methodology further. NiO, MnO and SrTiOs 
have local ground states which are well understood, and 
the spectra shown here have already been explained in 
the literature. New in the present paper is that the 
MLFT parameters (except for U and A) are not fitted to 
the experiment, but calculated ab initio. 

In the following subsections we first discuss x-ray ab- 
sorption (XAS) at the L2,s edge, which probes TM 2p 
to 3d excitations. Next we show TM 2p core- level x-ray 
photoemission (XPS) experiments on an impurity sys- 
tem. Both 2p XAS and 2p XPS excite the same core 
states and the difference is that in x-ray absorption the 
electron is excited into the local 3(i-shell whereas photo- 
emission excites the core electron into vacuum states. 
The resulting spectra are very different. We then con- 
tinue with core excitations measured with inelastic x-ray 
scattering (IXS). From a theoretical point of view, in- 
elastic x-ray scattering of core to valence excitations and 
x-ray absorption of core to valence excitations is very 
similar. The initial and final states probed are the same. 
The difference is that whereas x-ray absorption is mainly 
caused by dipole transitions, inelastic x-ray diffraction is 
caused by multipole transitions determined by the length 
of the transferred momentum. Finally in the last sub- 
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FIG. 7: (color online) Comparison of the experimental (thick 
red) and MLFT (thin blue) TM 2p (L2,3 edge) core-level x-ray 
absorption spectra for SrTiOs, MnO, and NiO. The experi- 
mental SrTiOa spectra are reproduced from Uehara et al..^^ 
the MnO spectra from Csiszar et al.,^^ and the NiO spectra 
from D. Alders et al.¥^ 



section, we show inelastic x-ray scattering of d-d excita- 
tions in NiO. These spectra are particularly instructive as 
they allow for a relative straight-forward understanding 
on how the different interactions contribute to each mul- 
tiplet excitation. For pedagogical reasons we provide a 
brief introduction to each of the experimental techniques. 
A more thoroughly description of these techniques can be 
found in textbooks, e.g. those by De Groot and Kotani,'^ 
Stohr^^ Schulke.^3 



A. ^2,3 edge x-ray absorption 

X-ray absorption spectroscopy (XAS) at the TM ^2,3 
edge is a technique whereby a TM 2p core electron is 
excited into the 3d valence shell. The excitation energy 
is in the x-ray range and varies from ~ 440 eV for Ti 
to ~ 855 eV for Ni. The excitations are dipole allowed, 
which make them so intense that spectra with very little 
noise can be obtained. The spectra split into two set of 



peaks, the L3 and L2 edges, due to spin-orbit coupling in 
the TM 2p core hole. This results in a 2pj^^/2 level (L3 
edge) lying |C2p above a 2pj^i/2 level (L2 edge). For 
core levels the relativistic spin-orbit coupling is strong 
and element dependent: C2p=3.78 eV for Ti, 6.85 eV for 
Mn, and 11.50 eV for Ni. Hence, in Fig. [7] we see the 
L3 and L2 edges at 640 and 650 eV for MnO, and at 852 
and 870 eV for NiO. 

For SrTi03, |C2p is of similar size as the 2p-3d mul- 
tiplet splitting, i.e. the L2 and L3 edges overlap. The 
splitting within an L2 or L3 edge is due to the combined 
interaction of covalent ligand-field effects and Coulomb 
interactions between the 3d electrons and between the 2p 
core hole and the 3d electrons. This leads to the relatively 
involved spectra with many features as seen in Fig.[7[ 
Even for SrTi03 where one might be tempted to relate 
the four intense peaks in the 2p XAS spectrum to exci- 
tations from the 2pj=^/2 or 2pj=i/2 orbitals into either 
the t2g or the Cg orbitals, the intensity ratios (4 : 2 be- 
tween excitations from respectively 2^^=3/2 and 2pj=i/2 
core holes, and 3 : 2 between excitations to respectively 
t2g and Cg states) do not follow this one electron pic- 
ture: Assigning the peaks at 458, 460, 463 and 465 eV 
to excitations of the form 2^^=3/2 ^2^, 2^^=3/2 e^, 
2pj=i/2 and 2pj=i/2 eg respectively, would 

yield the intensity ratios 12 : 8 : 6 : 4, which are clearly 
not observed. On the other hand, starting from a 2p^3d^ 
final-state configuration in a cubic crystal field, does yield 
the correct intensities,^^ plus several small peaks. Our ab 
initio results shown in blue in the figure, confirm this in- 
terpretation. Within the atomic 2p^3d^ excitonic picture 
the interpretation in terms of t2g and Cg excitations of the 
Ls and L2 edge for the four peaks might still be a rea- 
sonable starting point, but one should realize that there 
is a considerable mixing between states due to Goulomb 
interactions. 

The cluster eigenstates cannot be represented by single 
Slater determinants. For correlated TM compounds, the 
spectral line-shape is governed by multiplet effects lead- 
ing to involved spectral functions, not obviously related 
to the density of states .1^^^ The spectra are therefore 
generally used as fingerprints which contain unique fea- 
tures resembling the local ground-state properties. The 
energy of the final state is determined by local atomic- 
like physics. The intensity with which each state can be 
reached depends, via the optical selection rules, on the 
ground-state symmetry and the polarization of the light. 
This can lead to large spectral changes for small changes 
in the ground- st at e.^^^*^^*^ For example, a p electron can 
only be excited into an dx2_y2 orbital with x or y polar- 
ized light, but not with z polarized light. Changing the 
orbital occupation can therefore lead to a strong polar- 
ization dependence which for certain multiplets can be 
as strong as 100%.^^ Due to the strong TM 2p spin-orbit 
coupling the XAS spectr a are also sensitive to the spin 
of the ground-state.E^ 

Theoretically, as well as experimentally, one finds that 
the monopole part of the 2p-3d Coulomb interaction is 
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larger than that of the 3d-3d interaction.l^IE^ This 
leads to strongly bound excitons at the TM 1/2,3 edge and 
allows one to describe the spectra using MLFT. Besides 
these excitonic states, also excitations into non-bound 
states are possible. Such excitations essentially probe 
the conduction bands of the compound. For NiO both 
excitations are clearly visible in the experimental spectra 
shown in the bottom panel of Fig.[7| The excitonic bound 
states give rise to sharp excitations which extend upwards 
from 852 eV; they are seen to agree very well with our 
MLFT spectra. At 865 eV the experimental NiO spectra 
show an edge jump where the cross section for photon ab- 
sorption increases discontinuously. This is the onset of L3 
excitations into the conduction-band continuum without 
formation of bound-states. These continuum excitations 
are of course not reproduced with MLFT. The L2 excita- 
tions into bound excitons give rise to the sharp features 
starting at 870 eV and captured by our MLFT. Around 
856 eV there is a slight disagreement between the theo- 
retical and experimental Ni-I/3 edge spectra which might 
be due to the neglect of the Ni 45 orbitals in our cluster 
basis set. 

Looking at NiO, MnO and SrTi03, one may notice 
that our calculations reproduce the low-energy parts of 
the spectra better than the high-energy parts. The for- 
mer are most excitonic and therefore best described by 
the small basis set in the cluster. It may furthermore be 
noticed that not only the edge-jumps are absent in the 
calculation, but also the interference between the exci- 
tonic excitations of the L2 edge and the continuum states 
of the Ls edge. These interference effects give rise to 
Fano like line-shapes present in the experiment, but not 
in the theory. The effect is relative small as the interfer- 
ence between 2^^=3/2 and 2pj=i/2 states is forbidden in 
many channels. There is, however, a substantial mixing 
of core-states due to Coulomb interactions, which could 
be the main reason for the interference effects between 
continuum and excitonic states of the L3 and L2 edge. 

Nevertheless, the agreement between MLFT and ex- 
periment is rather good in Fig.[7| for all three TM ox- 
ides. This agreement is as good as -or even better 
than- that obtained for calculations in which all standard 
MLFT parameters are optimized to give the best fit to 
experiment. 9i 94 ^j^g^-g ^^^e many parameters in such a 
calculation and finding the best fit is not trivial. The use 
of ab initio values for an otherwise equivalent MLFT cal- 
culation can therefore be of great help to interpret x-ray 
absorption spectra and thus also in the understanding of 
elastic resonant x-ray diffraction (RXD) spectra, and 
the resonant energy dependence of resonant inelastic x- 
ray scattering (RIXS).^^^*^ For systems with lower local 
symmetry, the number of parameters is even larger, and 
so is the need for values determined ab initio. 

Compared with other ab initio methods used for the 
calculation of the L2.3 edge s of correlated transition- 
metal compound^^^J'^^'^^HI^ the current method pre- 
forms well. For d^ compounds, i.e. band in- 
sulators like SrTi03, very powerful methods based 
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FIG. 8: (color online) Comparison of the experimental (thick 
red) and MLFT (thin blue) Ni 2p core-level photoemission 
spectra of NiO in the impurity limit. The experimental spec- 
tra are reproduced from Altieri et a/.P^ 

on multiple-scatte ring formalisms j ^^*^^ ' or the Bethe- 
Salpeter equations,^^^ are available. For Mott-Hubbard 
or charge-transfer insulators, ab initio configuration- 
interaction calculations of high quality have been per- 
formed for finite-sized cluster s.'^^Mll q^j. MLFT method 
has the advantage that its one-electron basis functions 
exactly describe the relevant bands for the infinite crys- 
tal and at the same time localize so well that one can 
afford to include correlations beyond the LDA for merely 
the TM d orbitals. This allows for very efficient, but still 
accurate many-body calculations in the framework of the 
well-studied multiplet ligand field theory. 



B. 2p core-level photoemission of Nio.03Mgo.97O 

Core-level photoemission is uninteresting from a one- 
electron point of view. Core levels are atomic like, 
have no momentum-dependent dispersions, and there- 
fore delta-peaked densities of states. Accordingly, 2p 
photoemission is expected to yield two spin-orbit-split 
peaks with intensity ratio 2:1. However, the exper- 
imental spectra^^^ from Ni 2p core-level photoemission 
in Nio.03Mgo.97O exhibit four distinct spectral features, 
as shown in Fig.lsl Emission from the Ni 2^3/2 level 
gives rise to the structure between 852 eV and 868 eV 
and emission from the 2pi^2 level to the structure be- 
tween 870 eV and 886 eV. The structure between 860 
and 868 eV originates from multiplet excitations with 
main character 2p^^^3d^ (with the underbar indicating a 

hole) while the peak centered at 855 eV belongs to an ex- 
citation with main final-state character 2p 3d^L. It is 

— 3/2 

the strong Coulomb attraction between the 2p core hole 
and the 3d electron which causes the latter state to have 
lower energy than the former, and thus screens the core 
hole by driving charge in from the ligand. Between 870 
and 886 eV, this spectrum of screened and un-screened 
states is repeated, but now for excitations from the 2pi/2 
core level. Our MLFT spectra agree well with the exper- 
imental spectra^^^ and with MLFT calculations for Ni 
impurities in MgO with fitted parameters .^^^ The result- 



ing interpretations of the experiment are the same. 

These Ni 2p core-level photoemission spectra are strik- 
ingly different from the previously considered Ni 2p x-ray 
absorption spectra in the bottom panel of Fig.jTl^^ In 
core-level x-ray absorption, a TM core electron is excited 
into the 3(i shell of the same atom, whereby the sample 
remains locally neutral. The electron can bind with 
the core-hole left behind and thereby screen the core- 
hole potential. This gives rise to the strong excitonic 
peaks seen in the x-ray absorption spectrum. In core- 
level photoemission, a core electron is emitted from the 
sample (excited into the vacuum) and can therefore not 
screen the core hole left behind. The core- hole is either 
screened by the surrounding ligands or left un-screened, 
which gives rise to higher-energy excitations. 

Photoemission spectra are generally not excitonic. It 
might therefore seem strange to use MLFT to calcu- 
late those spectra. Nevertheless, experience has shown 
that many features of photoemission from correlated 
transition-metal compounds can be captured by full mul- 
tiplet theory for a local cluster. Photoemission combined 
with cluster calculations has contributed greatly to our 
understanding of correlated TM and RE compounds.^^ 
The influence of non-local screening, i.e. the effect of 
the material being a solid and not a single impurity, 
has been studied experimentally by comparing the core- 
level photoemission from TM impurities with that from 
the TM compounds The main features of the Ni 2p 
photoemission spectra from Nio.03Mgo.97O are the same 
as from NiO. The largest bulk effect is a splitting of the 
peaks at 855 and 873 eV. 

Important progress in understanding bulk valence 
photoemission from NiO has been made recently by 
solving the LDA O p Ni d Wannier-orbital Hub- 
bard model in the dy namical mean-field approximation 
(LDA+DMFT^ISIII03] 

and also by using the variational 
cluster approximation.fl^ How important correlations 
between different Ni sites is, and therefore how impor- 
tant the inclusion of dynamical non-local screening ef- 
fects are,1^ remains an open question. On a different 
level, MLFT is able to reproduce a substantial part of 
the photoemission spectra, even though these are not ex- 
citonic. 



C. Non-resonant inelastic x-ray scattering at the 
M2,3 edge 

In subsection A, we compared experiments and MLFT 
for core level x-ray absorption (XAS). Now we shall dis- 
cuss core-level spectra obtained with a technique which 
from a theoretical point of view is very similar to x-ray 
absorption, namely inelastic x-ray scattering (IXS). In 
XAS at the TM M2,3 edge, a TM 3p core state is ex- 
cited to a ?>d conduction state by absorption of a pho- 
ton. The same excitation can be made when a photon 
is scattered inelastically and only part of its energy is 
absorbed. 106 108 rpj^^ major difference between XAS 
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FIG. 9: (color online) Comparison of the experimental (thick 
red) and theoretical (thin blue) Mn 3p (M2,3 edge) non- 
resonant inelastic x-ray scattering spectra of MnO at large 
momentum transfer (q). The two panels give spectra for dif- 
ferent directions of momentum transfer. They exhibit the 
generalized natural linear dichroism present for an octupole 
transition in cubic symmetry. The experimental spectra are 
reproduced from Gordon et al.}^^^ 



and IXS is that, for the former, the energy of the photon 
has to equal the absorption edge, whereas for the latter, 
the energy of the photon should be (much) higher than 
the absorption edge since only a fraction of its energy 
is absorbed. For XAS at the TM M2,3 edge, the lead- 
ing interaction is of dipole character, i.e. one can use 
the long wave-length limit. For IXS, the transferred mo- 
mentum can be selected by looking at different scatter- 
ing angles and energies: for small momentum transfers, 
dipole transitions are measured and for larger momen- 
tum transfers, octupole transitions .^^^^ In Fig. [o] we show 
the non-resonant IXS at the M2,3 edge in MnO for high 
magnitude momentum transfers where octupole transi- 
tions are the strongest. 

There is a clear difference between the ^2,3 (2p to 
and M2,3 (Sp to 3(i) edges in MnO. One reason is that 
the spin-orbit coupling constant for 3p is much smaller 
than for 2p, e.g. for Mn, ("2^ = 6.85 and Csp = 0-77 
eV. The splitting between the M3 and M2 edges is thus 
much smaller (not resolved in the experimental spectra) 
than the splitting between the L3 and L2 edges. Another 
reason why the ^2,3 and the M2,3 edges look different 
comes from the fact that the 3p wave-function is larger 
than the 2p wave- function due to the extra node. This 
leads to a smaller monopole part of the Coulomb repul- 
sion and larger multipole interactions between the p and 
3<i orbitals. In general, for excitations within the same 
radial shell the multiplet splittings are larger than the 
excitonic binding energy. This has important conse- 
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FIG. 10: (color online) Comparison of the experimental (thick 
red) and MLFT (thin blue) non-resonant inelastic x-ray scat- 
tering intensity of low energy d-d excitations. The experimen- 
tal spectra are reproduced from Verbeni et al.r^ 



quences. The highest-energy multiplets of the M2^s exci- 
tations are pushed above the continuum edge and form 
broad resonances instead of sharp muhtplets. The low- 
energy multiplets, on the other-hand, are still sharp ex- 
citonic states. Due to the strict selection rules applicable 
to XAS and IXS, one can choose the experimental ge- 
ometry such that only particular excitations are allowed. 
The spectra shown in Fig. |9] are octupole dominated and 
only sensitive to the low-energy excitonic features in the 
spectra. One would not be able to reproduce the broad 
dipole resonances with MLFT. 

One of the beauties of octupole transitions is that they 
show dichroism in cubic symmetry.^^^ This can be seen 
in the two different panels of Fig.|9j and is well repro- 
duced by our theory. For a dipole transition one can 
not distinguish cubic from spherical symmetry. (A tran- 
sition of angular momentum L = 1 branches to a sin- 
gle irreducible representation (Ti^) in cubic symmetry.) 
An octupole transition, however, shows nice dichroism in 
cubic symmetry, whereby the momentum transfer q for 
IXS takes the place of the light polarization e in XAS. 
A transition of angular momentum L = 3 branches to 
three different irreducible representations in cubic sym- 
metry, namely Ti^, and A2u- As a consequence, the 
dichroic spectra can be used to obtain detailed informa- 
tion about the differences in bonding of t2g and Cg elec- 
trons. The shift in the spectral energy and the change in 
spectral weight for excitations with q either in the [111] 
(top panel) or [001] (bottom panel) direction is related 
to the different energy of the t2g and Cg electrons and the 
difference in occupation of these orbitals due to covalent 
bonding. 



D. Non-resonant inelastic x-ray scattering of d-d 
excitations 



This last section in which we compare MLFT with 
experiment deals with low-energy excitations without 
a core hole. It has recently been shown that surpris- 
ingly strong d-d or crystal-field excitations can be ob- 
served in NiO with non-resonant IXS for large momen- 



tum transfers liiittil^J These spectra contain similar in- 
formation as the weak d-d excitations inside the optical 
gap observed with optical spectroscopy.l^^^ The differ- 
ence between IXS and optics is that with optics these 
transitions, being even in parity, are forbidden and only 
become allowed by simultaneous excitation of a phonon 
and a crystal-field excitation. This makes a quantita- 
tive interpretation of optical d-d excitations involved. 
The interpretation of the non-resonant IXS is, on the 
other hand, str aight fo rward and allows for a quantita- 
tive comparison.^^SZEll 

we show the experimental^ and theoretical 



In Fig. 
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non-resonant IXS spectra for a powder of NiO at large 
momentum transfer (averaged over a transfer of 7.3 — 8.0 
A ~^). These spectra are governed by quadrupole and 
hexadecapole transitions between the 3d orbitals. The 
non-resonant IXS excitations are spin-conserving. Lo- 
cally the Ni ground-state configuration is d^ with the t2g 
orbitals fully occupied and the Cg orbitals half filled with 
(iS^) = 2, i.e. S = 1. In the one-electron picture, one 
can make a single excitation going from the t2g shell to 
the Cg shell, which has an experimental energy of about 
1.1 eV. This is the peak of T2g final state symmetry in 
the experiment. In principle one could also excite two 
t2g electrons simultaneously into the Cg sub-shell. This 
would give rise to a single peak at twice the energy. In a 
pure one-electron picture the double excitation is forbid- 
den because non-resonant IXS couples a single photon 
to a single electron. Using full multiplet theory, how- 
ever, both excitations have a finite intensity. This has 
to do with the strong t2g-eg multiplet interaction which 
mixes, for the excited states, the single Slater determi- 
nants. One even finds three peaks instead of two. The 
first peak indeed corresponds to an excitation of a single 
t2g electron into the Cg sub-shell. The second peak is 
roughly the simultaneous excitation of two t2g electrons 
into the Cg subshell. Finally, in order to understand the 
third peak, one should realize that the t2g {xy) electron 
is Coulomb repelled more from an Cg (x^ — y'^) electron 
than from an Cg (3z^ — l) electron because of the larger 
overlap of densities. This leads to multiplet splitting of 
the ^2^e^ states and to mixing of ^2^e^ and ^t^e^ states. 

One could also have understood the energy and num- 
ber of excitations by starting from spherical symmetry 
where Coulomb repulsion splits the 5* = 1 states into a 
lowest state of symmetry and an excited state of 
symmetry. In cubic symmetry the states branch into 
a ^A2 ground-state, a ^T2 first excited state, and a ^Ti 
second excited state. The state branches to a state 
of ^Ti symmetry, which can mix with the highest excited 
state branching from the ^F state. Such multiplet effects 
are hard to capture at the DFT level. Recent time de- 
pendent DFT calculations with the LDA+/7 functional 
do show Frenkel excitons {d-d excitations) within the op- 
tical gap, but they cannot reproduce the correct number 
of multiplet states.^^^ 

Let us finally have a closer look at the comparison be- 
tween the experimental and MLFT crystal-field excita- 
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tions in NiO. MLFT gets the lowest excitation (^Ts) 5% 
too low and the highest (^Ti) 10% too high. As the 
energy is mainly determined by one-electron inter- 
actions, we conclude that the Cg — t2g splitting due to 
covalency in our LDA based calculation is 5% underes- 
timated. At the same time, the multiplet splitting due 

(2) 

to the Coulomb repulsion, i.e. the values of the F^^^ and 

F^^ Slater integrals, are 10% overestimated. The later 
could be a result of neglecting the screening of the mul- 
tipole interactions, but not necessarily, because there are 
additional channels in which two 3d electrons can scatter 
into two higher excited states due to Coulomb repulsion. 
This gives rise to a multiplet-dependent screening, not 
easily described with a single screening parameter 



IV. CONCLUSIONS 

We have shown how multiplet ligand-field theory 
(MLFT) calculations can be based on ab-initio LDA 
solid-state calculations, in a sim ilar way as originally 
devised by Gunnarson et a/j ^^ l ^^ l and recently done for 
LDA+DMFT calculations. The resulting method could 
be named LDA+MLFT. The theory is very well suited 
for the calculation of local ground-state properties and 
excitonic spectra of correlated transition-metal and rare- 
earth compounds. Our TM d Wannier orbitals, which 
together with the O p Wannier orbitals span the LDA 
TM d- and O p-bands, are quite similar to atomic or- 
bitals, and this justifies many previous studies using the 
latter. 

We compared several experimental spectra (XAS, non- 
resonant IXS, PES) for SrTiOs, MnO, and NiO with our 
ab initio multiplet ligand-field theory and found overall 
satisfactory agreement, indicating that our ligand-field 
parameters are correct to better than 10%. The cova- 
lency seems to be slightly underestimated and the Slater 
integrals for the higher multipole interactions overesti- 
mated. The method is expected to provide insights to 
the local properties of transition-metal compounds with 
only modest computational efforts. 

We would like to thank Eva Pavarini, Ove Jepsen, 
and Olle Gunnarsson for fruitful discussions. Support by 
the Deutsche Forschergemeinschaft through FOR 1346 is 
gratefully acknowledged. 



Appendix A: Definition of O and TM orbitals, 
covalency and formal valence 

It is a general praxis to talk about Wi-d and 0-p or- 
bitals, even in a solid. The definition of local orbitals 
in a solid is not always clear. In this paper we choose 
such orbitals as generalized Wannier orbitals of a given 
local symmetry. For our materials, i.e. transition metal 
oxides, there are two useful definitions for the Wannier 
orbitals, depending on the size of the basis set used and 




FIG. 11: (color online) NiO band structure (left) and Wan- 
nier orbitals (right) for three different basis sets. Top panels: 
Including both the Ni d and O p orbitals. Middle panels: In- 
cluding only the Ni d orbitals. Bottom panels: Including only 
the O p orbitals. 



the energy bands they span. These two different defi- 
nitions are often used in an ambivalent way. Here we 
explain the definition of the two different set of orbitals, 
by the example of NiO. 
In Fig. 



11 we show in the top row the Ni d and O 



p Wannier orbitals which are most atomic like. Linear 
combinations of these 8 orbitals span the 8 bands shown 



on the left of the top row in Fig. 11 Due to covalency, the 
Ni d orbitals defined in this way are occupied by more 
than 8 electrons. At the same time, the O p Wannier 
orbitals are occupied by less than 6 electrons, i.e. they 
have holes. When forming band states, the Ni d and O 
p orbitals mix and part of the O p character ends above 
the Fermi energy. One should not think of these extra 
electrons or holes as mobile charge carriers. These partial 
occupations are just a result of the choice of the basis set 
used, which is different from the eigenbasis, the band 
states. 

A different choice of the Wannier orbitals, closer to 
the eigenfunctions, can be seen in the lower two rows of 
Fig. 11 Here we show from top to bottom the Ni d 
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(eg and t2g) and O p tiu orbitals. The Ni d orbitals are 
those 5 orbitals which together span the Ni d bands in 
the energy range from -3 to 2 eV. The O p orbitals are 
those 3 orbitals which together span the O p bands in 
the energy range from -8 to -3 eV. 

In order to distinguish the two different sets of Wannier 
orbitals we have throughout this paper used italic font for 
the more localized atomic like orbitals and roman font for 
the more extended orbitals. The difference in notation is 
quite subtle, but in almost all cases one can understand 
from the context which definition is meant. 

The two basis sets of either atomic like Ni d and O p or- 
bitals or more delocalized Ni d and O p orbitals span the 
same bands and can thus be expressed in terms of linear 
combinations of each other. The unitary transformation 
between the two sets of Wannier orbitals is such that it 
diagonalizes the covalent interaction between the Ni d 
and O p orbitals. The Ni d and O p orbitals interact, 
whereas the Ni d and O p orbitals are non interacting 
at the one particle or LDA level. The 5 Ni d orbitals 
span the 5 d bands exactly and the 3 O p orbitals spand 
the 3 O p bands exactly. The O p orbitals are bonding 
combinations of the TM d and O p orbitals. The TM d 
orbitals are anti bonding combinations of the TM d and 
O p orbitals. 

With the use of the Ni d (O p) Wannier orbitals, which 
span only the Ni d (O p) bands one can define the formal 
valence of Ni in NiO. It is common to state that O is 
2—, i.e. has an occupation of p^ and Ni is 2+, i.e. has 
an occupation of d^. If one counts the electrons in the 
Wannier orbitals that separately span the Ni or O bands 
one immediately reproduces the formal valence. The oc- 
cupation numbers are different if one looks at the more 
atomic like Ni d and O p orbitals. For these orbitals co- 
valence introduces holes in the O Wannier orbitals and 
extra electrons in the Ni Wannier orbitals. For the oxides 
described in the present paper, the occupations of the 
two kinds of Wannier orbitals are: NiO 3d^ ~ 2p^-^3d^-^^ 
MnO 3d^ ^ 2/•^3d^•^ and SrTiOa ^ {2p^-^) ^3d^-^ . 

Let us note that for our purposes, the Wannier d- 
orbitals are not sufficiently localized. Nevertheless, in 
early LDA+DMFT calculations which could handle only 
a few correlated orbitals, even more "downfolded" t2g 
or e^ Wannier-orbitals were used by necessity; they 
clearly exhibit the covalencies.l^ 



Appendix B: Computational details 

The selfconsistent LDAI^ L APW calculations were per- 
formed with the Wien2k code^^ using a plane- wave cut 
off of kmaxXRMT = 8, with Rmt the smallest MT- sphere 
radius and kmax the largest /c-vector. The A/^M TO calcu- 
lations were done with the Stuttgart cod4^^, had N=l^ 
and all p artial waves downfolded, except TM d and O 
^|47 h5Q | 52| rpj^^ LAPW warped potential (spherical inside 
the LAPW MT spheres) was least-squares fitted to an 
overlapping MT potential (OMT) with the recently de- 



veloped OMTA cod^ and was used in the A^MTO cal- 
culations. The radii of the hard screening spheres were 
70% of the OMT radii. 

The material-dependent settings are as follows: 

NiO. Space-group Fm-3m (225) a=4.177A, Ni at 
Wyckoff position 4a and O at 4b. MT radii for Ni 2.08ao, 
and for O 1.84ao. ao=0.5292 A, is the Bohr radius. 
OMT radii for Ni 2.2ao, for O 2.5ao, and for an addi- 
tional empty sphere at Wyckoff position 8c 1.6ao. The 
expansion energies were —5.2 and —1.2 eV. 

MnO. Space-group Fm-3m (225) a=4.4248A, Mn at 
Wyckoff position 4a, O at 4b. MT radii for Mn 2.20ao, 
and for O 1.95ao. OMT radii for Mn 2.3ao, for O 2.7ao, 
and for an additional empty sphere at Wyckoff position 
8c 1.7ao. Expansion energies —5.0 and —1.0 eV. 

SrTiOa. Space-group Pm-3m (221) a=3.905A, Ti at 
Wyckoff position la, Sr at position lb, and O at position 
3d. MT radii for Ti 2.32ao, for Sr 2.00ao, and for O 
1.36ao. OMT radii for Ti 2.4ao, for Sr 3.8ao and for O 
2.0ao. Expansion energies —2.6 and +1.5 eV. 

The A/'MTO band-structures and densities of states 
as presented in Fig.s [2] and [5] were calculated from the 
real-space TM 3d O 2p Wannier-orbital (tight-binding) 
representation of the LDA Hamiltonian, neglecting hops 
between sites more distant than 2.5a. The tight-binding 
parameters larger than 10 meV are presented in the lower 
panels of Table [T[ Only the nearest-neighbor hopping in- 
tegrals (bold faced) enter in the cluster calculations. In 
SrTiOs, the point symmetry of O is merely tetragonal so 
that the p orbital pointing towards the Ti atom is slightly 
different from those pointing perpendicular to the Ti-0 
bond (and e.g. towards Sr). For a discussion of the 
bonding between the O p and Sr d orbitals and how this 
changes the different O p Wannier functions see Pavarini 
et al.^^. Due to the two different types of O p Wannier 
functions, the relation between the hopping integrals and 
the cubic ligand-field parameters is slightly more involved 
than those valid for O in cubic symmetry and given in 
Sect. I. In general (for all symmetries), the ligand field 
parameters can be found by block tri-diagonalization of 
the tight-binding Hamiltonian of the cluster with respect 
to the TM d orbitals. For more details see Appendix [P] 

Wannier orbitals have tails on the neighboring sites, 
although most of the orbital weight is close to the nu- 
cleus at its center (Fig. 3|. The tails lead to long-ranged 
hopping integrals and their values are given in the lowest 
panels of Table [l[ It should be noticed that in order for 
MLFT to work properly with a basis set of LDA based 
Wannier orbitals, it is important to have TM 3d atomic- 
like character for r < 1.5 A, but it is not essential to have 
hopping limited to the first-nearest neighbors. 

The multipolar part of the Coulomb integrals is cal- 
culated by directly integrating the wave- functions. In 
order to obtain numerically stable integrals, the Wannier 
orbitals were expanded in radial wave-functions times 
spherical harmonics, an approximation for which a set 
of Slater integrals can be introduced. Sufficiently ac- 
curate results are obtained when different radial wave- 
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TABLE I: Upper panel: MLFT parameters obtained from LDA. The Slater integrals are obtained from the spherical averaged 
Wannier orbitals. Lower two panels: One-electron tight binding parameters as obtained from the LDA TM-d 0-p Wannier- 
orbital set. The hopping is from the first to the second orbital displaced by the vector [abc]. [000] denote on-site energies. 
Shown are only those values larger than 10 meV. The bold numbers enter in the MLFT calculations, the normal font longer 
range hopping integrals are truncated in the cluster approximation. For the p-p hopping in the [001] direction of SrTiOs the 
first two values listed concern hopping along an 0-Ti-O bond. The last three values concerns 0-0 hopping in the Sr-0 plane. 
The notation for the eg orbitals is such that (i^2 = d^z2_i and dxv = dx2_y2. All values in electron volt. 
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TABLE IL Multiplet ligand field theory parameters taken 
from experiment. Note that the experiments shown in this 
paper are not very sensitive to these parameters. For a more 
thoroughly discussion on these parameters we refer to the pa- 
pers by Bocquet (Ref. |78| or Tanaka (Ref. [79}. All values in 
electron volt. 



functions for the t2g and Cg orbitals are used. The core 
wave-functions, are calculated using the Hartree-Fock 
method. For reasons of space, Table I only gives the 
Slater integrals for the radial functions averaged over t2g 
and Cg] the difference between these integrals for NiO 
are given as an inset in Fig.[3| The spin-orbit coupling 
constants have been calculated using a spherical approx- 
imation, including only the local d character at the TM 
site. The resulting constant is the same for the t2g and 
Cg orbitals. 

The parameters fitted to experiments are shown in Ta- 
ble [nj These values are in good agreement with those in 
the literature.^l^ One should realize that since x-ray 
absorption involves a charge neutral excitation, it is not 
very sensitive to U and A. The experiments discussed 
in this paper were chosen to be most sensitive to the 
calculated values (Tab le . 

The 7VMT0 methodSEll constructs the basis set of Ni 
d plus O p localized orbitals by first constructing such a 
set for each of the N -\- 1 energies, eo,...,eAr, chosen to 



span the energy range of interest. In such a set of zeroth- 
order (A^=0) MTOs, each of the orbitals is a solution 
of Schrodinger's differential equation for the overlapping 
MT potential for the chosen energy, but has kinks at all 
Ni and O (but not e.g. Sr) hard spheres. Those hard 
spheres are chosen to be a bit smaller than touching and 
not to coincide with a node of the radial wave- function. 
The Ni dxy OMTO, for instance, is now defined by the 
hard-sphere boundary condition that all its p-projections 
on all O spheres and all its (i-projections on all Ni spheres, 
except dxy on the own Ni sphere, vanish. And equiv- 
alently for the other members of the OMTO basis set. 
This hard-sphere boundary condition is what localizes 
the OMTOs, unless there are wave-functions at the cho- 
sen energy with main characters different from those of 
the OMTOs in the basis. The condition that each OMTO 
solves Schrodinger's equation, except for kinks in the Ni 
d and O p channels, means that each OMTO is smooth 
in all other channels. This is accomplished by construct- 
ing that set of wave-equation solutions in the hard-sphere 
interstitial, the so-called screened spherical waves, whose 
phase shifts are the hard-sphere ones for all Ni d and O 
p channels, except the eigen-channel, and has the proper 
phase shifts for all other channels (such as Ni s and Sr 
d) . The screened spherical waves then gets augmented in- 
side the overlapping MT spheres to become the OMTOs. 
Finally, the + 1 different OMTO basis sets are con- 
tracted into one, the A^MTO set, which spans the solu- 
tions of Schrodinger's equation at all A^ + 1 energies. The 
A/'MTOs have discontinuities in merely the (2A^ + l)th 
derivatives at the hard spheres and are therefore smooth 
if > 0. The contraction (A'-ization) delocalizes the 
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A/'MTO to a degree which depends on how much the 
neighboring OMTOs vary over the + 1 energies. This 
is so, because for an energy-independent set of orbit als, 
the energy dependence of a radial Schrodinger-equation 
solution must be provided by the tails of the neighboring 
orbitals.E^The delocalization is further enhanced by sym- 
metrical orthonormalization of the A^MTOs into Wan- 
nier orbitals, and this depends on the overlap between 
neighboring A'MTOs. Nevertheless, as seen in Fig.|2j our 
Ni d plus O p Wannier orbitals are as localized as can 
be expected, and -in fact- much better than Wannier 
orbitals derived from a large set of energy- mdependent 
orbitals.E^ 



Appendix C: Double-counting correction 

DFT in the LDA already contains a large part of the 
local Coulomb interactions. These interactions are in- 
cluded once more in the MLFT calculations. For MLFT 
based on the LDA potential and Wannier orbitals one 
should therefore take care not to double count such inter- 
actions. We differentiate rigidly between the monopole 
and the multipole parts of all Coulomb interactions. This 
idea is based on the experimental observation that the 
monopole part of the Coulomb interaction (U) is largely 
screened, from ~ 25 eV to ^ 7.3 eV in NiO for example. 
The multiplet splitting, determined by the multipole part 
of the Coulomb interaction is, however, only slightly re- 
duced from the splitting one expects based on atomic 
values. This has for example been observed in Auger 
spectroscopy for the elemental 3d metals.^ 

The monopole part of the Coulomb interaction ([/) as 
well as the spherical part of the on-site energy (A) we 
fit to the experiment and double-counting for the spher- 
ical part is therefore not an issue. In order to prevent 
double counting of the multipole part of the Coulomb 
repulsion, the LDA calculations are done with a warped 
LDA potential; i.e. within the MT sphere only the spher- 
ical part of the potential is included. In order to check 
how this influences the LDA band-structure we compare 
in Fig. 12 the band-structures of NiO (top) and SrTiOs 
(bottom) calculated with the full LDA potential (thick) 
and warped MT potential (thin). Concentrating first on 
the NiO bands, we see that both calculations agree within 
basically the line- width for all bands and k- vectors, ex- 
cept for the t2g bands. Those bands are shifted down- 
wards in the warped-MT calculations by a momentum 
independent value of about 220 meV. This effect nicely 
illustrates the problem of double counting. The orbital 
occupation of the Ni 3d orbitals within LDA is such that 
the t2g orbitals are fully occupied and the Cg orbitals are 
half filled. The local charge density is thus cubic. The 
Coulomb repulsion between two t2g orbitals is on average 
larger than between a t2g and eg orbital. This effect is 
well included in the LDA functional and related to the 
fact that the overlap of for example the density of the dxy 
and dxz orbitals is larger than the overlap of the density 




FIG. 12: (color online) Comparison of the LAPW band- 
structures and densities of states for the full LDA potential 
(thin lines) and the warped MT potential (thick lines) . When 
only a single line can be seen, the two band-structures over- 
lap within the line- width of the plot. The colors indicate the 
partial- wave characters inside the MT spheres. 



of the dxy and dsz'^_i orbitals. In MLFT calculations, 
such interactions are included in the Slater integrals. A 
MLFT calculation based on the full LDA potential would 
thus double count the multipolar interaction between the 
fully occupied t2g and the half filled eg shell. 

One option would be to include the full potential 
within the LDA calculations and then subtract the non- 
spherical part of the Coulomb repulsion, as included in 
the LDA functional, between the Wannier functions for 
which a full multiplet interaction is included in the MLFT 
calculations. In that case one should carefully analyze 
the occupation of each Wannier function in order to de- 
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termine the potential that has to be subtracted. We 
opted to not include the non-spherical interactions in 
the first place. This does mean that one also neglects 
the non-spherical part of the non-local Madelung poten- 
tial in the self-consistent LDA calculations. In order to 
correct for this, we calculate this potential from the self- 
consistent LDA charge density and added it afterwards. 
We found that the Coulomb potential which is double 
counted generally exceeds the non-spherical non-local po- 
tential by an order of magnitude. Our choice of not in- 
cluding the non-spherical interactions in the first place 
and treating the non-local non-spherical interactions as 
a correction after self-consistency has been reached thus 
provides an accurate self-consistent solution to the poten- 
tial needed in MLFT. Nevertheless, we expect that doing 
a full-potential calculation and subsequently subtracting 
the non-spherical part of the local Coulomb interaction 
according to the LDA functional will give very similar 
results. 

In the bottom panel of Fig.[l2] we compare the band- 
structures of SrTiOa calculated for the warped and full 
potentials. First of all, there is no clear shift of the t2g 
bands with respect to the Cg bands, presumably because 
Ti atom has a d^ configuration. But there are changes 
in the O p-derived bands. In NiO both the Ni and the O 
atoms have cubic point symmetry, but in SrTiOs the O 
environment is tetragonal as was mentioned above. For 
SrTiOs one can see a larger difference between the full 
and warped-potential calculations than for NiO. The in- 
terpretation for SrTiOs is less straight forward because 
several effects come together. Due to Ti-0 covalency 
there is some non-spherical, predominantly e^^-derived 
charge on the Ti atom from the O 2p band. The lo- 
cal non-spherical potential due to this charge should not 
be included when doing MLFT calculations. There is 
however also a non-cubic potential at the O site that 
shifts the O bands. This potential should be included 
when doing the MLFT calculations. Note that the latter 
potential is not included in the self-consistent warped- 
potential calculations, but we add it later, before doing 
the MLFT calculations. 



Appendix D: Basis size reduction and the creation 
of Ligand Orbitals by blocktridiagonalization 

Embedded cluster calculations, like MLFT or DMFT, 
contain a few correlated orbitals coupled to a large set of 
uncorrelated orbitals. Within such calculations the size 
of the Hilbert space can be reduced enormously by creat- 
ing appropriate linear combinations of the uncorrelated 
orbitals. Within our examples these are the O p orbitals, 
which are combined to Ligand orbitals. Without intro- 
ducing ligand orbitals for a TMOe cluster, there would 
be (18 + 5) X 2 = 46 spin orbitals in the one-electron 
basis. With a filling of 8 + 6 x 6 = 44 electrons for the 
NiOe d^-cluster, this results in 46!/(44! x 2!) 10^ states 
in the many-electron basis. For the TiOe d^-cluster rep- 



resenting SrTiOs, the filling would be 6 x 6 =36 and 
thus result in 46!/(36! x 10!) - 4 x 10^ states in the 
many-electron basis. The introduction of L orbitals, how- 
ever, reduces the number of one-electron basis functions 
to 20, whereby the many-electron Hilbert space reduces 
to 20!/(18! X 2!) = 190 for a d^ and to 20!/(10! x 10!) - 
2 X 10^ states for a d^ configuration. This reduction in 
the number of many-electron basis functions by factors of 
respectively ||| ^ ^ W\ (46lio)! ^ leads 

to a crucial gain of computational convenience. Either 
basis set can be used to calculate ground-state properties 
and spectral functions because the matrices are sparse. 
But diagonalization of a matrix with dimension 4 x 10^ 
requires large computational resources whereas diagonal- 
ization and evaluation of spectral functions of a sparse 
matrix with dimension 2 x 10^ can be done using stan- 
dard libraries on modern desktop computers. One may 
obtain a further reduction in the number of stored basis 
states by including only those which are important for 
representing the actual wave- function, (see Appendix [E| 
for details). 

Ligand orbitals are normally obtained by symmetry 
considerations.'^ The rotation properties of the TM d 
orbitals should be the same as the linear combination of 
the O p orbitals with which this orbital makes a covalent 
bond. These symmetry considerations can be extended 
to a simple mathematical procedure, valid in all point 
group symmetries. From DFT the one particle Hamil- 
tonian for an extended cluster is known on a basis of 
the central TM d orbitals and the neighbor O p orbitals. 
Using a block Lanczos routine one can create a unitary 
transformation of the p orbitals such that the one parti- 
cle Hamiltonian has a blocktridiagonal form. The basis 
of the central TM d orbitals is not changed. In cubic 
symmetry the tridiagonalization results in a transformed 
Hamiltonian whereby each TM d orbital couples to one 
Ligand orbital. For lower symmetries each d orbital cou- 
ples to maximal 5 Ligand orbitals. The ligand orbitals 
can couple to an other set of ligand orbitals, ad infini- 
tum. Covalence in a tridiagonal representation tends to 
converge fast, justifying the inclusion of only a single Lig- 
and shell in MLFT. 

The introduction of Ligand orbitals is not restricted 
to fully occupied shells, like the O p shell. For example 
for SrTiOs one could include besides the O p Ligand or- 
bitals also the Sr d Ligand orbitals. Care has to be taken 
how for such a system the Ligand orbitals are defined. If 
one creates a single Ligand shell for both the O p and 
Sr d orbitals by blocktridiagonalization as described in 
the previous paragraph one would obtain Ligand orbitals 
that are always partially occupied. This results in a very 
large many particle basis set and is unpractical. In such 
a case it is better to first diagonalize the non-interacting 
Hamiltonian describing the interactions between the Lig- 
and orbitals. Based on the onsite energies one then cre- 
ates two different Ligand shells, one for the occupied, or 
valence orbitals and one for the unoccupied or conduction 
orbitals.EH 
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Note that a similar procedure can be used for DMFT 
calculations using a Lanczos impurity solver. Doing so 
enhances the calculation speed and allows one to increase 
the number of bath sites (number of discretization sites 
used to represent the Anderson impurity model used in 
DMFT) leading to much more continues spectral func- 
tions. 



Appendix E: Exact diagonalization and Lanczos 
algorithm 

The MLFT ground-state and s pectra l calculations are 
done using a Lanczos algorithm !^^ * ^^ ! The calculations 
start with a random vector {ipo) in the basis of the d'^ 
and d^^^L configurations, whereby n is the number of 
d electrons (0 for Ti, 5 for Mn and 8 for Ni) and L 
represents a single hole in the Ligand shell. Although 
this starting point is slightly worse than the DFT sin- 
gle Slater determinant ground-state, which prescribes a 
specific mixture of d and L states, it does contain the 
correct symmetry states. Thereby convergence is so fast 
that the starting point really does not matter that much. 
Given a negative definite Hamiltonian, the wave-function 
V^i = HipQ has a larger overlap with the ground-state 
wave function than the wave- function ?/^o- By repeat- 
edly acting with the Hamiltonian on the random starting 



function and normalizing the wave-function in-between 
(V^n+i = Hi/jn/lHi/jnl), one converges to the ground- 
state. This procedure can be speeded up considerably by 
creating a tri-diagonal matrix of the Hamiltonian in the 
basis of i/jn^ with the additional constrained that V^n+i is 
orthogonal to tpn- The tri-diagonal matrix in the so called 
Krylov basis can be diagonalized with the use of dense 
matrix methods. Having found the ground-state within 
the basis of the d^ plus d^^^L configurations, we remove 
the basis functions not needed to represent the ground- 
state wave function from the basis and extend the basis 
set by acting with the Hamiltonian on the wave- function. 
This creates basis states belonging to the d^^'^I? con- 
figuration. Within this new basis set the ground-state is 
found and the procedure of truncating and extending the 
basis set is repeated. The whole process is repeated un- 
til convergence is reached. Excited states are calculated 
by repeatedly orthogonalizing the wave-function to the 
eigenstates already found. The algorithm as described 
here allows one to always keep relatively small basis sets. 

Spectral functions are calculated by acting with the 
transition operator on the ground-state wave- function. 
The resulting function is then used as a starting vector 
for the creation of a tri-diagonal matrix in a Krylov ba- 
sis. The spectral function of a tri-diagonal matrix can be 
expressed in terms of a continued fraction. 
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